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摘要 : 为 了 深入 分 析 和 探讨 昆虫 发 育 与 环境 温度 的 关系 , 在 恒温 (15 ~37%C ) 和 交替 变温 (12/18 ~34/40C ) 下 测定 了 
棉铃 虫 Helicoverpa armigera 肾 的 发 育 历 期 (d) ,分 别 用 线性 模型 和 非 线性 模型 ( Logan 模型 、 Lactin 模型 和 王 氏 模型 ) 
拟 合 其 发 育 率 (1/d) 数 据 。 结 果 表 明 ,这 3 个 非 线性 模型 能 更 准确 地 描述 发 育 率 与 温度 之 间 的 曲线 关系 ,判定 系数 
(R°) Æ 0.9878 ~0.9991 之 间 。 对 全 部 观测 数据 的 进一步 研究 表明 ,只 要 有 6 个 分 布 合适 的 观测 数据 ,就 可 以 用 这 些 
非 线性 模型 获得 相当 满意 的 佑 测 效果 。 Me UU 用 非 线 性 模型 预测 的 昆虫 发 育 率 可 能 失真 。 分 
析 了 晴 在 恒温 和 变温 下 发 育 率 差异 的 可 能 原因 ,讨论 了 应 用 这 3 个 非 线性 模型 预测 旺 期 发 育 的 优点 和 缺点 ,指出 用 
非 线 性 模型 取代 线性 日 . 度 模型 进 和 Fa te NO RS Aie tS Am RUA TETE c 
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Estimating developmental rates of Helicoverpa armigera ( Lepidoptera: 


Noctuidae) pupae at constant and alternating temperatures by nonlinear models 
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Abstract: Developmental durations (d) of Helicoverpa armigera ( Hübner) pupae were determined at 
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constant temperatures from 15 to 37°C and alternating temperatures from 12/18 to 34/40'C for 
comprehensive analyses on insect development in relation to temperature. Both linear model and nonlinear 
models ( Logan model, Lactin model and Wang model) were used to fit the observed data, and these 
nonlinear models gave a more accurate description of temperature-dependent developmental rates with the 
coefficients of determination ( R^) ranging from 0. 9878 to 0. 9991. Further fitting various numbers of the 
observed data showed that so long as six cases of them distributed in the proper way, these nonlinear models 
could estimate pupal developmental rate with high satisfaction. On the other hand, lack of the data at high 
temperature would give a false prediction. A possible cause for difference between developmental rates of 
the pupae at constant and alternating temperatures is illustrated. The strong and weak points of the three 
nonlinear models when applied to estimate the pupal developmental rate are commented. The results suggest 
that the nonlinear models may be of significant value in improving prediction of insect pest occurrence as 
well as rearing and management of beneficial insects as alternative to linear degree-day model. 
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温度 是 影响 昆虫 生长 发 育 的 主要 环境 因素 之 


加 快 ,然后 加 速 增加 到 最 大 值 ,此 后 随 痢 温度 的 继续 


Raven 的 快慢 通常 用 其 发 育 历 期 (d) 或 
历 期 的 倒数 (1/d)， 即 每 天 完成 的 发 育 率 
( developmental rate) 表示 。 许 多 研究 表明 , 昆虫 发 
育 率 与 温度 之 间 是 一 种 复杂 的 曲线 关系 :从 发 育 起 
点 温度 开始 ,昆虫 的 发 育 率 随 着 温度 的 升 高 而 逐渐 
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公式 是 早期 最 常见 的 两 种 方法 ,它们 虽然 应 用 简便 ， 
但 并 不 能 全 面 反映 昆虫 在 其 能 存活 的 整个 温 区 内 的 
发 育 率 变化 情况 。20 世纪 70 年 代 非 线性 回归 模型 
的 研发 (Logan et al., 1976) 是 该 领域 研究 的 一 个 重 
要 突破 ,此 后 有 多 个 非 线性 模型 相继 问世 ( Sharpe 
and DeMichele ,1977 ; Taylor, 1981 ; 王 如 松 等 ，1982 ; 
Hilbert and Logan, 1983; Lamb et al.,1984;Lactin et 
al.,1995 ; Briere et al., 1999) 。 这 些 模 型 都 能 在 不 同 
程度 上 更 全 面 地 描述 昆虫 或 螨 类 发 育 率 与 温度 之 间 
曲线 关系 。 

为 了 深入 分 析 和 探讨 昆虫 发 育 与 环境 温度 的 关 
系 ,我 们 利用 棉铃 虫 Helicoverpa armigera ( Hübner) 师 
在 系列 恒温 和 变温 下 的 发 育 率 测定 数据 ,选用 常见 的 3 
个 非 线 性 模型 分 别 进行 拟 合 , 均 获 得 相当 满意 的 结 采 。 
在 此 基础 上 ,进一步 研究 了 用 非 线性 模型 预测 昆虫 发 
育 率 所 需要 的 最 少 观测 数据 例 数 。 根 据 拟 合 结果 , 讨 
论 和 评价 了 这 3 个 非 线性 模型 的 优 、 缺 点 ,指出 了 用 非 
线性 模型 符 代 线性 日 度 模 型 进行 害虫 发 生 预 测 和 益 
虫 人 饲养 管理 的 合理 性 和 必要 性 。 


1 材料 和 方法 


1.1 HER 

以 室内 连续 饲养 多 代 的 棉铃 虫 肾 作为 研究 对 
象 。 选 择 棉铃 虫 肾 是 因为 在 用 人 工 饲 料 饲养 的 例 行 
实验 (27Y ,L14 : D10) 中 经 和 常 有 大 量 多 余 的 老龄 幼 
忠 , 稍 作 处 理 即 可 利用 。 同 时 , 肾 期 不 取 食 ,无 需 频 
繁 检 查 , 化 肾 和 羽化 时 间 也 易于 确定 。 
1.2 温度 处 理 

试验 在 人 工 气 候 箱 内 进行 ,温度 误差 +1C ,分 
别 在 恒温 和 交替 变温 下 测定 肾 的 发 育 历 期 。 变 温 幅 
度 6% ,每 天 在 比 相 应 恒温 高 4% 的 条 件 下 处 理 8 h, 
然后 移 至 比 相 应 恒温 低 27C 的 条 件 下 处 理 16 ph, 使 
日 均 温 度 等 于 相应 的 恒温 。 

将 老 训 幼虫 置 于 含 十 (含水 量 约 8% ) 的 指 形 管 
内 , 待 其 入土 后 即 转移 到 设 定 的 温度 下 处 理 ,每 天 
9 :00 检查 化 肾 或 羽化 数 一 次 并 区 分 性 别 。 处 理 期 
间 无 光照 ,及 时 在 管内 补充 水 分 ,保持 适当 湿度 。 每 
处 理 供 试 的 两 性 晴 数 均 在 50 头 左 右 。 整 个 实验 的 
工作 量 很 大 ,历时 数 年 完成 。 
1.3 数据 的 统计 处 理 和 模型 的 拟 合 

肾 在 不 同 恒温 和 变温 下 的 发 育 历 期 均 以 能 羽化 
的 成 虫 (不 论 正常 与 否 ) 数 为 基础 进行 统计 分 析 , 由 
于 上 肉 、 雄 肾 的 历 期 差异 显著 ,分 别 统 计 , 取 历 期 的 倒 


数 (1/d) 作 为 每 天 完成 的 发 育 率 。 通 过 SPSS 10.0 
软件 绘制 温度 -发 育 率 关系 的 散 点 图 ,根据 散 点 图 不 
同 区 有 段 的 特征 和 特定 观测 值 的 位 点 设 定 非 线 性 模型 
中 有 关 参 数 的 起 始 值 ,选用 以 下 3 个 在 文献 中 引用 
较 多 、 应 用 效果 较 好 的 非 线 性 模型 拟 合 观测 数据 。 

Logan 模型 :Logan 等 (1976) 在 分 析 温 度 对 迈 
IHi Tetranychus mcdanieli 发 育 的 有 影 啊 时 ,将 两 者 之 
间 关 系 的 曲线 分 为 两 部 分 ,一 部 分 是 从 发 育 起 点 温 
度 以 上 某 个 低温 下 的 发 育 率 随 温 度 升 高 而 增加 到 最 
大 值 时 的 单调 上 升 型 指数 曲线 ,为 一 部 分 是 发 育 率 
达到 最 大 值 后 随 温度 继续 升 高 而 急剧 降低 的 “边缘 
层 ” 曲 线 , 分 别 赋予 两 者 数学 函数 表达 式 , 然 后 整合 
成 单一 的 非 线 性 模型 , 它 在 文献 中 的 常见 形式 是 : 

r(T) = wiexp(pT) - exp[pT,, - (Ta - T)/AT |] 

(1) 
式 中 :r (T) 表示 温度 7 了 时 的 发 育 率 ,y 表示 发 育 起 
点 温度 以 上 某 个 基础 温度 下 的 发 育 率 ,p 表示 发 育 
率 随 温度 变化 的 指数 增长 率 ,T, 表示 致死 高 温 ,AT 
表示 DAJE” 宽度 。Logan 模型 是 文献 中 引用 最 多 
的 非 线 性 模型 之 一 。 

Lactin 模型 :Lactin 等 (1995 ) 仔细 审视 了 Logan 
模型 后 ,对 该 模型 进行 了 两 处 修改 ,一 是 删除 多 余 的 
参数 由 ,二 是 加 入 一 个 新 的 参数 “入 ”: 

r(T) = exp(pT) - exp| pT, - (T, - T)/AT | +à 

(2) 
式 中 :参数 入 表示 估 测 的 发 育 起 点 温度 。 近 年 来 ， 
Lactin 模型 已 被 广泛 用 于 分 析 多 种 昆虫 发 育 率 与 温 
度 的 关系 (Briere et al. , 1999; Roy et al., 2002; 
Fantinou et al., 2003 ; Gould et al., 2005; Davis et al., 
2006; Chen et al., 2006; Tsoukanas et al., 2006; 
Kalaitzaki et al., 2007 ; Zamani et al., 2007) , 

王 氏 模型 :这 是 我 国学 者 王 如 松 等 (1982 ) 在 归 
纳 前 人 有 关 理 论 和 经 验 公 式 的 基础 上 推导 出 的 非 线 
性 模型 : 


WT) = {1 -expl - (T - 


K 

1 + expl - r(T - T,) | 
T,)/8]]1 i1 - expl - CT, - T)/8]] (3) 
式 中 :V(7T) ,r, T, 和 6 分 别 相当 于 Logan 模型 中 的 > 
(T) ,p,T, 和 AT7; 天 = 最 大 发 育 率 (内 ) 的 2 倍 ,To = 
最 大 发 育 率 ( WV) 时 的 温度 ,TL = 发 育 起 点 温度 。 我 
国学 者 用 该 模型 分 析 了 温度 对 两 种 蚜虫 ( 刘 树 生 和 
dz 989) 和 小 菜 蛾 (Liu et al., 2002) 发 育 的 影 
啊 , 均 获得 满意 的 预测 效果 。 
1.4 评价 指标 

用 上 述 模型 通过 迭代 法 对 观测 数据 进行 拟 合 ， 
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根据 校正 判定 系数 (adusted coefficient of 
determination, R^), 残 差 平方 和 (residual sums of 
squares, RSS) 以 及 预测 的 发 育 起 点 温度 、 最 适 温度 
和 上 限 温度 比较 和 评价 拟 合 效 果 。 


2 结果 和 分 析 


2.1 温度 对 棉铃 虫 肾 期 发 育 的 影响 
棉铃 虫 肾 期 发 育 对 温度 变化 的 反应 很 敏感 ,在 
恒温 31~ 32C 或 变温 29/35%C 以 下 ,发 育 历 期 随 温 


表 1 棉铃 虫 肾 在 恒温 和 


度 的 升 高 而 缩短 ,在 高 于 这 些 温度 时 , 发育 历 期 则 随 
着 温度 的 升 高 而 趋 于 延长 ,尤其 是 高 于 恒温 34% 或 
变温 31/37%C 时 ,发 育 历 期 大 幅度 延长 。 尽 管 进行 
了 重复 试验 ,但 在 恒温 32 ~ 34% 或 变温 29/35 
~ 32/38'C 之 间 仍 然 未 能 找 出 发 育 率 由 升 高 转 为 降 
低 的 明显 转折 点 ( 表 1 ,图 1)。 在 恒温 38Y 或 变温 
35/41*C 下 , 肾 虽 然 能 存活 不 同时 间 ,但 最 终 无 一 羽 
化 ,全 部 死亡 。 





变温 下 的 发 育 历 期 (d) 


Table 1 Developmental duration in day of pupae of the cotton bollworm, Helicoverpa 


armigera ( Hiibner) at constant and alternating temperatures 


BUE CC) 2 3 
Temperature CT AT CT AT 

14 43.48 +0.29 (49) 48. 78 +0.18 (50) 
15 41.24 +0.20 (50) a 37.32 +0.25 (49) b 45.45 +0.15 (41) a 41.68 +0.24 (45) b 
16 37.16 +0.25 (32) a 32.89 «0.17 (35) b 41.15 +0.15 (31) a 35.09 +0.22 (35) b 
17 32.47 +0.21 (32) a 28.09 «0.19 (33) b 36.76 +0.21 (32) a 31.25 +0.19 (33) b 
18 26.22 +0.17 (27) a 23.76 +0.16 (32) b 30.12 +0.23 (30) a 26.53 «0.12 (32) b 
19 22.71 +0.20 (35) a 19.27 +0.19 (33) b 25.85 +0.12 (40) a 21.55 «0.16 (33) b 
20 19.01 +0.14 (35) a 17.44 +0.17 (29) b 21.75 +0.10 (35) a 19.69 «0.12 (32) b 
21 16.79 «0.10 (36) a 15.56 «0.10 (27) b 19.27 +0.10 (35) a 17.89 «0.09 (40) b 
22 15.63 +0.09 (32) a 14.50 +0.15 (33) b 17.76 x0.14 (32) a 15.83 «0.18 (33) b 
23 14.27 +0.15 (34) a 12.58 «0.10 (29) b 16.18 +0.12 (29) a 14.16 +0.14 (30) b 
24 12.48 +0.09 (29) a 11.43 +0.10 (35) b 14.01 +0.08 (37) a 12.82 «0.12 (36) b 
25 11.30 +0.11 (39) a 10.48 +0.08 (33) b 12.66 «0.09 (34) a 11.48 +0.11 (33) b 
26 10.56 «0.07 (39) a 9.93 +0.08 (30) b 11.56 «0.04 (39) a 10.89 «0.07 (38) b 
27 9.69 «0.09 (50) a 9.27 +0.08 (43) b 10.67 +0.08 (35) a 10.28 «0.09 (44) b 
28 9.14 +0.06 (33) a 8.56 +0.04 (42) b 10.03 +0.11 (31) a 9.60 +0.07 (28) b 
29 8.24 +0.06 (33) a 8.27 +0.07 (39) a 9.31 +0.08 (33) a 9.21 +0.08 (39) a 
30 8.02 +0.09 (27) a 7.94 +0.08 (34) a 8.90 «0.08 (30) a 8.76 +0.07 (37) b 
31 7.43 +0.09 (35) a 7.72 +0.05 (41) b 8.14 +0.07 (37) a 8.51 «0.10 (34) b 
32 7.27 +0.11 (32) a 7.79 «0.05 (32) b 7.89 +0.05 (26) a 8.53 +0.10 (31) b 
33 7.35 +0.05 (31) a 7.81 +0.11 (34) b 7.95 +0.06 (29) a 8.51 +0.11 (31) b 
34 7.49 +0.13 (32) a 7.97 «0.09 (32) b 8.03 +0.09 (32) a 8.71 +0.06 (36) b 
35 7.63 +0.09 (52) a 8.37 +0.08 (40) b 8.16 +0.06 (49) a 8.84 +0.08 (36) b 
36 7.90 +0.08 (74) a 8.85 «0.08 (40) b 9.10 «0.06 (98) a 9.82 «0.12 (42) b 
37 10.03 «0.22 (17) 12.35 «0.37 (11) 


表 中 数据 为 平均 值 土 标准 误 。Data in the table are mean € SE. CT; 恒温 Constant temperature; AT; 变温 ,在 比 相应 恒温 高 4% 下 处 理 8 h, 然后 
在 降低 6% 的 条 件 下 处 理 16 h, 使 日 均 温 等 于 相应 的 恒温 Alternating temperature, exposed to 4% higher than the corresponding CT for 8 h, then to 
6*C lower than it for 16 h, so that the daily average temperature is equal to the respective CT. 括号 内 数字 是 羽化 的 成 虫 数 , 数 据 后 有 不 同 字 母 表示 
差异 显著 ( P< 0.05, :检验 )。The numerals in parentheses are number of adults emerged, and the data followed by different letters indicate 


significant difference at P< 0.05 by t test. 下 同 The same for the following tables and figures. 
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图 1 是 以 表 1 中 发 育 历 期 的 倒数 (1/d), 即 每 
天 完成 的 发 育 率 与 温度 关系 的 散 操 图, 它 更 直观 地 
描述 了 温度 对 肾 期 发 育 率 影 响 的 全 貌 :同样 温度 下 ， 
雌 晴 发 育 率 普 志 大 于 雄师 ;在 平均 温度 相同 时 , 师 在 


发 育 率 Developmental rate (1/d) 
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恒温 和 变温 下 的 发 育 率 的 差异 很 大 ,只 有 在 恒温 29 

~30C (9 ) 和 变温 27/33%C ( 8 ) 下 , 晴 发 育 率 几乎 
相等 ; 低 于 此 温度 时 变温 下 的 发 育 率 均 显 著 大 于 和 恒 
温 下 的 相应 值 ,高 于 此 温度 时 则 相反 。 
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平均 温度 Mean temperature (°C) 
图 1 棉 伶 虫 晴 发 育 率 与 温度 关系 的 散 点 图 


Fig. 1 Scatter diagram of developmental rates of Helicoverpa armigera pupae in relation to temperature 


2.2 棉铃 虫 晴 期 发 育 率 与 温度 关系 的 拟 合 

2.2.1 用 线性 模型 和 logistic 曲线 模型 拟 合 :棉铃 
虫 晴 在 恒温 和 变温 下 的 发 育 率 变化 ,只 在 一 定 温 度 
区 间 与 环境 温度 呈 近 似 直 线 关 系 (图 1) 。 在 整个 实 
验 温 度 范围 内 ,用 线性 模型 和 logistic 曲线 模型 拟 合 
两 者 之 间 的 关系 都 不 能 令 人 满意 。 以 雌 晴 在 恒温 下 
的 发 育 率 变化 为 例 ,用 线性 模型 和 logistic 曲线 模型 
拟 合 时 的 R? 分 别 约 为 0.87 和 0.78 ( 表 2) ,都 不 能 
反映 发 育 率 在 高 温 区 随 温度 升 高 而 降低 的 特征 (图 
2), 


o 观测 值 Observed data 
~ 直线 Linear 


Logistic 








育 
Developmental rate (1/d) 


Y Ti 13 15 17 19 21 23 25 27 31 33 35 37 
温度 Temperature (C) 
图 2 用 线性 模型 和 logistic EAIA HETE 
雌 晴 在 恒温 下 的 发 育 率 
Fig. 2 Development rates of Helicoverpa armigera pupae ( 9 ) 


at constant temperature estimated linear and logistic models 


2.2.2 ”用 非 线性 回归 模型 拟 合 :用 Logan 模型 、 
Lactin 模型 和 王 氏 模型 拟 合 师 在 恒温 和 变温 下 发 育 
率 的 效果 都 相当 好 (图 3) ; R 均 达 到 或 大 于 0.99, 
按 尼 和 RSS 判断 ( 表 3) ,以 王 氏 模型 拟 合 效果 最 


佳 ;Logan 模型 不 能 预测 发 育 起 点 温度 (LT) ,用 王 氏 
模型 预测 的 LT 普遍 低 于 用 Lactin 模型 预测 的 LT; 
除 用 王 氏 模型 预测 上 肉 肾 在 变温 下 (DFA ) 的 致死 高 
温 (UT) 偏 高 外 ,3 个 模型 预测 的 最 适 发 育 温 度 
(OT) ftl UT 都 接近 于 观测 值 。 从 表 3 可 以 看 出 ,R 
与 RSS 的 关系 是 前 者 越 大 ,后 者 越 小 。 因 此 , 只 需 
R^ 即 可 判断 用 非 线性 模型 拟 合 结果 的 准确 性 。 
表 2 用 线性 模型 和 Logistic 曲线 模型 估 测 棉铃 虫 
肾 在 恒温 和 变温 下 的 发 育 率 
Table 2 Developmental rates of Helicoverpa armigera 
pupae at constant and alternative temperatures 


estimated by linear and logistic models 


模型 Model 预测 方程 Prediction equation R LT 


DFC 线性 Linear y 20. 005427x —0. 05272 0.86531 9-10 


Logistic y 21/(7.1429 + 929.76 x0. 8011* ) 0. 7766 - 


DMC 线性 Linear y =0. 004918x —0. 04846 0.8272 9-10 


Logistic y=1/(7.6923 +874. 39 x0.8098*) 0.7514 - 


DFA 线性 Linear y =0. 005307x — 0. 04658 0.9057 6/12 -7/13 


Logistic y =1/(7. 6923 +1 824. 39 x 0. 7631* ) 0. 7554 - 


DMA 线性 Linear y =0. 004892x — 0. 04433 0.9156 7/13 - 8/14 


Logistic y =1/(8. 3333 +1 330. 35 x 0. 7828" ) 0. 8648 - 


DFC, DMC: 分 别 表示 上 肉 、 雄 肾 在 恒温 下 的 发 育 率 Indicating 
developmental rates of female and male pupae at constant temperatures, 
respectively; DFA, DMA: 分 别 表示 上 肉 、 雄 肾 在 变温 下 的 发 育 率 
Indicating developmental rates of female and male pupae at alternating 
temperatures, respectively. LT; 发 育 起 点 温度 Threshold temperature 
(^C). 下 同 The same for the following tables. 
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图 3 JH Logan, Lactin 和 王 氏 模型 估 测 棉铃 虫 肾 在 恒温 和 变温 下 的 发 育 率 


Fig. 3 Developmental rates of Helicoverpa armigera pupae at constant and alternating temperatures estimated 
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10 13 16 19 22 25 28 31 34 37 40 8/14 
恒温 Constant temperature (°C) 


by Logan, Lactin and Wang models, respectively 
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Table 3 Accuracy of predicting developmental rates of Helicoverpa armigera pupae by the three nonlinear 


models fitted to observed data at different temperatures 


DFC DFA 
指标 Criteria 
Logan Lactin ER Wang Logan Lactin ER Wang 
R? 0. 9942 0. 9946 0. 9972 0. 9878 0. 9956 0. 9991 
RSS ( x10 ^?) 2. 0022 1. 8499 0. 9565 3. 8367 1. 3697 0. 2859 
LT 一 12 - 13 10 -11 一 9/15 - 10/16 2/8 -3/9 
OT 33 34 33 31/37 31/37 30/36 
UT 39 -40 38 -39 39 -40 36/42 -37/43 37/43 -38/44 40/46 -41/47 
DMC DMA 
指标 Criteria 
Logan Lactin ER Wang Logan Lactin ER Wang 
R? 0. 9965 0. 9949 0. 9985 0. 9903 0. 9961 0. 9991 
RSS (x10-4) 1. 0383 1. 5014 0.4370 2.5680 1.0262 0. 2259 
LT - 12 - 13 7-8 - 9/15 - 10/16 2/4 -3/3 
OT 33 34 33 31/37 31/37 31/37 
UT 38 -39 38 -39 38 -39 37/43 —38/44 37/43 -38/44 37/43 -38/44 


OT: 最 适 温度 Optimal temperature (0C) ; UT: 发 育 上 限 温 度 Upper threshold temperature (^C ). 下 同 The same for the following tables. 


2.3 观测 数据 例 数 对 用 非 线 性 模型 拟 合 结果 的 


影响 


表 1 列 出 了 晴 在 恒温 和 变温 下 发 育 历 期 的 观测 
数据 各 23 例 ,实验 温度 几乎 履 关 了 贤能 存活 的 整个 
温 区 ,获得 这 些 数据 既 费 时 又 费力 , 显然 并 非 上 策 。 


然而 ,这 些 资 料 给 我 们 提供 了 一 个 可 徘 的 平台 ,可 以 
从 中 任意 选取 大 干 例 数据 ,探讨 其 对 用 非 线 性 模型 
拟 合 结果 的 影 啊 。 

2.3.1 用 非 线 性 模型 拟 合 晴 期 发 育 率 所 需要 的 最 
少 观 测 数 据 例 数 : 按 均匀 分 布 原则 从 表 1 中 选取 不 
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ul 光 维 性 模型 要 | nz O 观测 值 Observed data 
同 数 据 例 数 ? 用 非 线 竹 模型 拟 H WARR 明 ? 只 要 0.18 M 被 拟 合 的 观察 值 Observed data to be fitted 
选取 的 数据 分 布 恰当 ,即使 只 有 6 例 数据 ,也 能 得 到 0.16 Tn 


— ER Wang 





相当 满意 的 预测 效果 ( 表 4) 。 以 雌 师 在 恒温 下 的 发 


育 率 (RFC ) 为 例 , 选取 在 15,20,25,29, 33 和 , R7 

37%C 下 测定 的 6 例 发 育 率 数据 ,用 Logan 模型 、 Ei. 

Lactin 模型 和 王 氏 模型 拟 合 得 到 的 R 均 接近 或 大 $ 006 

T 0.99 ,对 全 部 观测 数据 的 校正 判定 系数 ( Ru ) 在 ^ aM 

0.98 以 上 ,预测 的 LT,OT 和 UT( 表 4, 图 4) 也 与 用 Ln 

这 3 个 模型 拟 合 全 部 观测 数据 的 结果 (图 3, #3) s onm moa m 32 35 moa 

相似 。 图 4 用 Logan 模型 ，Lactin 模型 和 王 氏 模型 拟 合 6 例 
当选 取 的 观测 数据 减少 到 5 例 时 ,用 Logan 模型 观测 数据 估 测 棉铃 虫 肉 师 在 恒温 下 的 发 育 率 

拟 合 的 效果 依然 较 好 ,但 用 Lactin 模型 拟 合 晴 在 恒温 Fig. 4 Developmental rates of Helicoverpa armigera pupae 

下 发 育 率 (RFC,RMC ) 的 结果 失真 , 按 全 部 观测 数据 ( 9 ) at constant temperatures estimated by Logan, 

计算 ,其 预测 值 的 RA [4.0.90 £A. MERAIH Lactin and Wang models fitted to 6 cases of 

合 时 因数 据 太 少 ,迭代 程序 不 能 运行 ( 表 5)。 observed data respectively fitted to the six observed data 


表 4 用 3 个 非 线 性 模型 拟 合 6 例 棉 铃 虫 旺 发 育 率 数据 时 预测 结果 的 可 靠 性 
Table 4 Reliability of prediction by 3 nonlinear models fitted to six cases of temperature-dependent 


developmental rates in pupal stage of Helicoverpa armigera 


DFC* DMC* 
指标 Criteria —— 
Logan Lactin ER Wang Logan Lactin ER Wang 
R? 0. 9948 0. 9962 0. 9998 0. 9980 0. 9968 1. 0000 
Ra? 0.9919 0.9900 0.9919 0.9935 0. 9899 0. 9962 
LT - 11 -12 8-9 - 11 -12 9 -10 
OT 33 34 33 33 34 34 
UT 39 — 40 38 —39 39 —40 38 —39 38 —39 38 -39 
PR T T T T T T 
DFA” DMA” 
指标 Criteria 一 一 
Logan Lactin ER Wang Logan Lactin ER Wang 
R 0. 9862 0.9970 0.9998 0. 9898 0.9979 0.9999 
Ra? 0. 9828 0.9911 0.9967 0.9874 0.9928 0.9980 
LT - 8/14 -9/15 3/9 4/10 - 8/14 -9/15 5/11 -6/12 
OT 31/37 31/37 30/36 31/37 31/37 30/36 
UT 37/43 -38/44 37/43 -38/44 40/46 - 41/47 37/43 ~ 38/44 37/43 —38/44 39/45 — 40/46 
PR T T T T T T 


Rn”: 对 全 部 观测 值 的 校正 判定 系数 Adjusted coefficient of determination for all observed data ( Kvalseth, 1985). PR: 对 全 部 观测 数据 预测 的 可 靠 
性 Predictive reliability for the all the observed data; T = IE True; F = AE False. 下 同 The same for the following tables. a: 拟 合 的 温度 -发 育 率 
数据 Fitted to the cases at 15, 20, 25, 29, 33 and 377€ ; b; 拟 合 的 温度 -发 育 率 数据 Fitted to the cases at 12/18, 17/23, 22/28 , 26/32 , 30/36 and 
34/40. 


2.3.2 ”观测 数据 分 布 对 拟 合 结 采 的 影响 : 按 2.3.1 当选 取 的 数据 偏 于 低温 区 、 测 定 的 最 高 温度 为 
节 的 方法 选取 6 例 温度 -发 育 率 数据 ,但 这 些 数 据 的 。 327023 29/35'C , 即 观测 的 发 育 率 最 大 时 ,用 3 个 非 
分 布 俩 于 低温 或 高 温 一 侧 时 ,用 非 线性 模型 拟 合 的 线性 模型 对 选 定 的 6 例 数据 的 拟 合 效果 仍然 较 好 
结果 差别 很 大 。 (R? »0.97) ,但 对 全 部 数据 的 拟 合 结果 , 除 个 别 例 
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外 ,普遍 不 佳 (Ra < 0.93) ,甚至 出 现 Ra < 0 的 
极度 失真 现象 ( 表 6) S OS EN YE THIS P ASI SII HU 
合 结果 清 楚 地 表明 ,失真 现象 主要 出 现在 缺乏 选 人 
数据 的 32C 以 上 高 温 区 (图 5)。 

当选 取 的 数据 偏 于 高 温 区 、 测 定 的 最 低温 度 为 
恒温 20C 或 变温 17/23%C 时 ,用 3 个 非 线 性 模型 对 
选 定 的 6 例 数据 和 全 部 观测 数据 的 拟 合 结果 都 相当 
满意 ,其 R 二 0.99,R 三 0.98 ,但 用 王 氏 模型 预测 


HJ LT 显然 太 低 ( 表 7) 。 
O ”观测 值 Observed data 
0.20 e UG Observed data to be fitted 
Ul] o -T---- Logan P 
—-- -—-- - Lactin PP 
0.16 \ 


d 
| 一 
N 


0.08 





14 17 20 23 26 29 32 35 38 
温度 Temperature (°C) 


图 5 用 Logan 模型 ，Lactin 模型 和 王 氏 模型 拟 合 指定 的 
6 例 发 育 率 观测 值 预测 棉铃 虫 肉 肾 在 恒温 下 的 发 育 率 
Fig. 5 Development rates of Helicoverpa armigera pupae 
( 9 ) at constant temperature predicted by Logan, Lactin 

and Wang models fitted to 6 indicated cases 


of observed data, respectively 
3 讨论 


3.1 发 育 速率 与 发 育 率 
研究 昆虫 发 育 与 温度 的 关系 时 普遍 采用 昆虫 发 
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育 历 期 的 倒数 (1Xd) 描述 其 发 育 的 快慢 ,后 者 的 英 
文 是 “developmental rate" , 其 中 文 一 直 沿 用 “发 育 速 
率 ” 表 述 。 仔 细 推 项 起 来 ,这 种 表述 方式 值得 商 榨 : 
“速率 ”是 指 “ 运 动 的 物体 在 单位 时 间 内 所 经 过 的 
距离 ”现代 汉语 词典 ,pl 086) ,其 中 的 时 间 和 距离 
都 应 该 有 具体 度量 单位 。 发 育 历 期 的 倒数 (1/d) 并 
不 具备 这 种 属性 , 它 是 指 每 天 完成 发 育 的 比率 。 
此 ,本 文 将 这 个 术语 简称 为 “发 育 率 ”。 
3.2 ”变温 和 恒温 对 昆虫 发 育 的 不 同 影响 

昆虫 在 自然 条 件 下 所 处 的 环境 温度 是 不 断 变化 
的 ,而 研究 温度 对 昆虫 发 育 影响 的 试验 一 般 都 在 恒 
温 下 进行 。 许 多 人 研究 表明 ,在 恒温 下 的 试验 结果 并 
不 能 完全 反映 昆虫 在 自然 条 件 下 的 真实 情况 
( Back, 1983 ; Hanula et al., 1987 ; xy Bj ^E fI EZ, 
1989 ; Liu et al., 1995 ,2002; Fantinou et al., 2003; 
Davis et al., 2006) 。 我 们 的 测定 结果 也 清楚 地 表明 
恒温 和 变温 对 棉铃 虫 肾 发 育 的 不 同 影响 :在 低温 区 ， 
肾 在 变温 下 的 发 育 比 在 恒温 下 快 ,在 高 温 区 则 相反 ， 
只 有 在 个 别 中 间 温 度 下 , 晴 在 恒温 和 变温 下 的 发 育 
率 才 比较 接近 (图 1) 。 寻 致 这 种 现象 的 原因 可 能 是 
多 方面 的 ,昆虫 在 恒温 下 发 育 率 与 温度 之 间 的 曲线 
关系 显然 是 其 中 的 一 个 重要 因素 。 以 表 1 中 的 数据 
为 例 ,简要 谢 析 如 下 :在 恒温 24 ~37% 之 间 ,实测 的 
晴 发 育 率 与 温度 之 间 大 体 上 呈 抛 物 线 关系 (图 6)， 
按照 试验 设计 ,平均 温度 为 313 的 变温 是 在 29 
350C 下 分 别处 理 16 h 和 8 h。 这 样 的 变温 处 理 下 ， 
肾 每 天 实际 完成 的 累计 发 育 率 (0. 1311 224 x8 + 
0.1214 +24 x16 = 0.1246) 显然 比 在 31% 恒温 下 


RS 用 3 个 非 线性 模型 拟 合 5 例 棉 铃 虫 晴 发 育 率 数据 时 预测 结果 的 可 靠 性 


Table 5 Reliability of prediction by the three nonlinear models fitted to five cases of temperature-dependent 


developmental rates in pupal stage of Helicoverpa armigera 


DFC? DMC* 
指 标 Criteria 

Logan Lactin Logan Lactin 
R? 0.9975 0. 9988 0. 9989 0. 9967 
Raj 0. 9881 0. 9026 0. 9947 0. 8989 
LT 一 11 - 12 一 11 - 12 

OT 33 36 33 36 
UT 39 -40 37 -38 38 -39 37 -38 

PR T F T F 


DFAP DMA” 
Logan Lactin Logan Lactin 
0. 9906 0. 9988 0. 9983 0. 9983 
0. 9807 0. 9872 0. 9859 0. 9903 
- 8/14 -9/15 - 9/15 - 10/16 
31/37 31/37 31/37 31/37 
37/43 - 36/42 — 36/42 — 35/41 — 
38/44 37/43 37/43 36/42 
T T T T 


a: 拟 合 的 温度 -发 育 率 数 据 Fitted to the cases at 15, 20, 26, 32 and 37€ ; b; 拟 合 的 温度 -发 育 率 数据 Fitted to the cases at 12/18, 17/23 , 23/29, 
29/35 and 34/40'C. 用 王 氏 模型 拟 合 5 例 数 据 时 ， 由 于 数据 太 少 ,迭代 不 能 运行 。There are fewer cases than the iteration to operate when Wang 


model is applied to the 3 cases of data. 
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表 6 用 3 个 非 线 性 模型 拟 合 棉 铃 虫 晴 在 6 个 指定 温度 下 发 育 率 数据 时 预测 结果 的 可 靠 性 
Table 6 Reliability of predicting developmental rates in pupal stage of Helicoverpa armigera by the 


three nonlinear models fitted to the six selected observed data 


DFC? DMC? 
指标 Criteria 
Logan Lactin ER Wang Logan Lactin ER Wang 
R? 0. 9706 0.9979 0.9997 0.9982 0.9951 0.9998 
Ra? 0. 0857 0. 7034 0.7881 0. 9867 0. 6375 0. 6360 
LT - 11 -12 2 -3 - 12 - 13 « -10 
OT 43 38 > 50 33 43 > 50 
UT 44 —45 39 -40 > 50 38 —39 45 -46 > 50 
PR F F F T F F 
DFAP DMA” 
指标 Criteria 
Logan Lactin ER Wang Logan Lactin ER Wang 
R? 0. 9964 0. 9968 0. 9997 0.9975 0. 9964 0. 9998 
Ra? 0. 4701 < 0 0. 9076 0.7405 < 0 0. 9280 
LT - 8/14 -9/15 « -8/ -14 - 9/15 — 10/16 < -8& -14 
OT 29/35 29/35 » 48/54 29/35 29/35 » 48/54 
UT 34/40 -35/41 29/35 -30/36 » 48/54 34/40 —35/41 29/35 — 30/36 » 48/54 
PR F F F F F F 


a: 拟 合 的 温度 -发 育 率 数据 Fitted to the cases at 15, 18, 21, 24, 28 and 32% ; b: 拟 合 的 温度 -发 育 率 数 据 Fitted to the cases at 12/18, 15/21, 
18/24, 21/27 , 25/31 and 29/35'C. 


表 7 用 3 个 非 线 性 模型 拟 合 棉铃 虫 肾 在 6 个 
Table 7 Reliability of predicting developmental rates in pupal stage of Helicoverpa armigera by the 


= 28 
A XE im 


three nonlinear models fitted to the six selected observed data 


度 下 发 育 率 数据 时 预测 结果 的 可 靠 性 


DFC* DMC* 
指标 Criteria 
Logan Lactin ER Wang Logan Lactin ER Wang 
R 0. 9967 0.9937 0.9985 0. 9988 0. 9966 0. 9988 
Ra? 0. 9907 0. 9884 0.9915 0.9932 0. 9888 0.9926 
LT - 13 - 14 « -30 - 13 ~14 « -30 
OT 33 33 33 33 34 33 
UT 39 — 40 39 —40 40 - 41 38 —39 38 —39 38 -39 
PR T T T T T T 
DFAP DMA” 
指标 Criteria 
Logan Lactin ER Wang Logan Lactin ER Wang 
R? 0. 9870 0. 9938 0. 9980 0. 9903 0. 9970 0. 9987 
Ra? 0. 9785 0.9889 0.9970 0.9812 0.9901 0.9955 
LT - 10/16 -11/17 < -28/ -34 - 10/16 - 11/17 « -28/ -34 
OT 31/37 31/37 30/36 31/37 31/37 30/36 
UT 37/43 —38/44 38/44 —39/45 42/48 -43/49 37/43 —38/44 37/43 —38/44 40/46 —41/47 
PR T T T T T T 


a; 拟 合 的 观测 数据 Fitted to the cases at 20, 23, 26, 29, 33 and 377€ ; b; 拟 合 的 观测 数据 Fitted to the cases at 17/23, 20/26, 23/29, 26/34, 


30/36 and 34/40*C. 
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温度 Temperature (°C) 
图 6 雌 晴 在 恒温 (CT) 和 变温 (AT) 下 发 育 率 差异 原因 的 图 解 


Fig. 6 Diagram for differences between developmental rates of Helicoverpa armigera pupae ( 9 ) 


at constant ( CT) and alternating ( AT) temperatures 
DAT : 29/35'C 变温 下 的 计算 值 Developmental rate calculated at alternating temperature (29/35 'C ). 


的 观测 值 小 得 多 ( 见 图 6 中 箭头 所 示 ) RZ, ME 
低温 区 变温 下 的 发 育 必 然 比 在 恒温 下 更 快 些 。 
此 ,根据 恒温 下 的 实验 数据 分 析 和 预测 昆虫 在 目 然 
变温 条 件 下 的 发 育 进度 时 ,要 根据 具体 情况 进行 适 
当 修 正 。 至 于 变温 本 身 是 否 卫 接 有 影响 昆虫 的 发 育 ， 
已 有 学 者 探讨 (Wagner et al., 1984 ; XI E RU i ^£ 
多 ,1989 ;Liu et al., 1995) ,不 再 歼 述 。 
3.3 温度 对 昆虫 发 育 影 响 的 预测 

在 昆虫 能 存活 和 发 育 的 温度 范围 内 ,其 发 育 率 
与 温度 之 间 是 一 种 复杂 的 曲线 关系 (图 1) ,用 简单 
的 线性 模型 或 logistic 模型 来 描述 这 种 复杂 的 关系 
无 疑 会 产生 相当 大 的 偏差 (图 2) ,与 之 相 比 ,用 非 线 
性 模型 描述 和 预测 温度 对 昆虫 发 育 的 影响 要 准确 得 
多 (图 3, 表 3)。 可 徘 的 实验 数据 是 应 用 数学 模型 
进行 分 析 的 基础 。 从 理论 上 说 ,这 样 的 数据 越 多 ,用 
非 线 性 模型 预测 的 结果 越 可 靠 , 但 在 实践 中 ,大 家 都 
希望 用 较 少 的 时 间 获 得 足够 的 必要 数据 。 本 文 从 全 
部 观测 数据 中 分 别 选 取 硅 干 例 ,用 非 线 性 模型 拟 合 
的 结果 表明 ,只 要 实验 设计 得 当 , 获 得 6 组 分 布 合适 
的 温度 -发 育 率 数据 ,用 这 3 个 非 线 性 模型 拟 合 都 可 
以 得 到 满意 的 预测 效果 (图 4, 表 4,6) ,又 好 又 快 地 
完成 科研 任务 。 限 于 篇 幅 , 本 文 只 能 列举 选取 厂 干 
例 观 测 数据 ,用 非 线 性 模型 拟 合 的 结果 。 实 际 情况 
复杂 得 多 ,有 时 用 某 个 非 线性 模型 拟 合 测定 数据 可 
失效 ,在 这 种 情况 下 ,需要 检查 或 适当 改变 模型 中 
些 参数 的 初始 值 ,也 可 以 改 用 其 他 非 线性 模型 拟 
,只 要 测定 数据 分 布 合理 ,一 般 都 能 获得 满意 的 


n 


Anm BS 
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R^ 是 评判 用 非 线性 模型 拟 合 实验 数据 优 劣 的 
主要 依据 ,这 只 是 指 对 既定 数据 的 拟 合 情 况 。 在 试 
验 数 据 有 限 、 分 布 不 尽 合理 时 ,只 根据 R 来 判断 


拟 合 和 预测 效果 未 必 可 靠 。 如 表 6 所 示 , 用 3 个 非 
线性 模型 拟 合 6 组 数据 时 ,表面 看 来 效果 很 好 ,R 
都 在 0.97 以 上 ,但 由 于 这 些 数 据 的 分 布 偏 于 低温 一 
侧 , 根 据 全 部 观测 值 和 各 自 对 应 预测 值 计 算 的 判定 
系数 Ru , 除 少数 例外 ,普遍 不 佳 。 为 了 全 面 评价 在 
数据 较 少 时 用 非 线性 模型 预测 棉铃 虫 肾 的 发 育 情 
况 , 本 文采 用 了 两 个 辅助 指标 :R 和 发 育 起 点 温度 
(LT) 、 最 适 温度 (OT) 、 上 限 温 度 (UT)( 表 4,7)。 
在 通常 情况 下 ,无 法 计算 Ru ,只 能 根据 模型 输出 的 
LT,OT 和 UT, 辅 之 以 预测 图 形 ( 图 4,5) ,判断 其 预 
测 结果 是 否 可 靠 。 

3.4 3 个 非 线 性 模型 应 用 效果 的 评价 

从 模型 的 复杂 性 来 看 , 王 氏 模型 居 首 ( 含 6 个 
参数 ) ,其 余 2 个 模型 只 含 4 个 参数 。 用 非 线性 模型 
进行 迭代 运算 前 ,必须 恰当 地 设置 有 关 模 型 中 参数 
的 初始 值 , 这 就 意味 着 模型 中 涉及 的 参数 越 多 ,对 观 
测 值 的 数量 及 其 分 布 情况 的 要 求 越 严格 。 例 如 ,用 
Logan 模型 拟 合 5 例 观 测 数据 也 能 获得 较 好 结果 ， 
在 同样 条 件 下 , 王 氏 模型 则 不 能 运行 ( 表 5)。 

f R^ 的 大 小 判断 , 王 氏 模 型 优 于 其 他 模型 。 然 
而 , 表 3 所 列 的 12 个 R^ ffi XE 0.9878 ~0.9991 之 
H, 差异 并 不 大 ,预测 效果 都 很 好 。 

LT 是 我 们 希望 预测 的 种 群 参数 之 一 。 由 于 昆 
虫 在 临近 LT 的 低温 下 发 育 极 慢 , 常 导 致 其 最 终 死 
亡 ,实际 上 不 可 能 通过 实验 准确 测定 昆虫 的 LT 
( Herrera et al., 2005) , 只 能 利用 数学 模型 估计 其 
LT。 如 表 3 所 示 ,Logan 模型 不 能 估 测 LT; 与 Lactin 
模型 相 比 ,用 王 氏 模型 预测 的 LT 更 低 些 ,该 模型 估 
测 肉 、 雄 肾 在 变温 下 的 LT 分 别 为 2/8 ~ 3/9Y 和 
-2/4 ~ -3/3 ,考虑 到 棉铃 虫 在 晴 期 的 清 育 特性 ， 
根据 经 验 和 常识 来 判断 ,这 种 估计 值 很 难 接 受 。 在 
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表 7 中 ,用 王 氏 模型 预测 晴 发 育 的 Rao 都 大 于 
0. 99 ,但 预测 的 LT 竟然 低 于 -30% 或 -28/ -349%C ， 
更 是 难以 置信 。 

综 上 所 述 , 这 3 个 非 线 性 模型 各 有 长 短 : Logan 
模型 和 Lactin 模型 比较 简单 ,应 用 方便 ,可 以 获得 满 
意 的 预测 结果 ,前 者 的 缺点 是 不 能 佑 计 LT, 后 者 舍 
计 的 LT 比较 接近 现实 。 王 氏 模 型 比较 复杂 , 预测 
结果 更 精确 些 , 但 有 时 预测 的 LT 过 低 , 很 难 接受 。 

近年 来 许多 研究 表明 ,与 目前 广泛 应 用 的 线性 
模型 相 比 ,用 非 线性 模型 预测 昆虫 的 发 育 更 准确 、 
更 可 靠 ( 刘 树 生 和 对 学 多 ，1989; Macro et al., 
1997; Lopez et al., 2001; Liu et al., 2002; Fantinou 
et al., 2003; Herrera et al., 2005) 。 现 在 的 研究 结 
末 表 明 :只 要 实验 设计 得 当 , 有 6 个 分 布 合适 的 观测 
数据 , 即 可 选用 适当 的 非 线 性 模型 预测 昆虫 的 发 育 
( 表 4,5; 图 4) ,配合 线性 模型 通常 也 要 5 个 观测 数 
据 ; 就 计算 过 程 而 言 ,用 非 线 性 模型 预测 昆虫 的 发 育 
似乎 更 复杂 些 , 但 计算 机 的 普及 已 使 繁琐 的 计算 过 
程 变 得 相当 便捷 。 人 至 于 SPSS 统计 分 析 软 件 的 应 用 
和 操作 ,已 有 几 个 版 本 的 中 文教 材 可 供 参 考 。 我 们 
要 切实 贯彻 科学 发 展 观 的 精神 ,及 时 更 新 观念 ,逐步 
采用 非 线性 模型 分 析 和 预测 温度 对 昆虫 发 育 的 影 
啊 ,使 我 国 的 害虫 预测 预报 和 益 虫 的 饲养 管理 工作 
提高 到 一 个 新 的 水 平 。 
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